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Abstract 

We study the dnamics of small magnetic cantilevers, either made from Si covered by a mag- 
netic film or entirely ferromagnetic ones. The magnetomechanical torques are found to cause line 
splittings in ferromagnetic resonance spectra and magnetization reversal facilitated by mechanical 
degrees of freedom. We show that the magnetomechanical torques can extend the limits of detect- 
ing and exciting motion at the nanoscale. A "nanomotor" described here effectively transforms 
rf magnetic fields into mechanical oscillations. We furthermore propose to integrate mechanical 
oscillators into magnetoelectronic devices that make use of current-induced spin-transfer torques. 
This opens new possibilities for electric transducers of nanomechanical motion. 

PACS numbers: 
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I. INTRODUCTION 



Micro- and nanoelectromechanical systems (NEMS)^ allow nanoscale control over as 
small sensors with spatial resolution on an atomic scale 3 operating at frequencies in the 
GHz range.— NEMS detection of extremely small forces corresponding to biomolecular 
interactions^ and mass changes corresponding to single molecules^ have been reported. Ef- 
forts have been made to merge the field of nanomechanics with that of nanoscale mag- 
netism. The success in the detection of a single electron spin by magnetic resonance is a 
good example.— A nuclear spin sensor that is based on imposing a coherent motion on nuclear 
spins has been proposed.^ Here we focus on the possibilities provided by the coupled motion 
of the strain field of a cantilever and the magnetization of a ferromagnet.^i Magnetization 
reversal in small magnetic clusters^ is usually realized by external magnetic field a 12 ' 1 ^' 14 ' 1 ^ or 
polarized spin currents employing spin-transfer torques . 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 Along alternative 
mechanisms to switch the magnetization such as employing antiferromagnets^ and time- 
dependent magnetic fields 2 ^ we suggest the effect of mechanically assisted magnetization 
reversal.^ The fastest mechanical device reported in the literature is operated with the help 
of magnetomotive forces at GHz frequencies.— We propose to actuate such systems by de- 
magnetizing currents that provide coupling between the mechanical and the magnetization 
motion most efficiently in the GHz range. Coupling of the magnetization motion to an elec- 
tric circuit via spin-transfer torques and magnetoresistance then opens the possibility of high 
frequency transducers of nanomechanical motion. In this paper, we report consequences for 
the magnetovibrational coupling such as splitting of the ferromagnetic resonance (FMR) 
spectrum and magnetization reversal facilitated by the mechanical degrees of freedom. We 
argue that a predominantly ferromagnetic cantilever provides stronger coupling and, when 
integrated into a magnetoelectronic circuit, offers new functionalities for both mechanical 
and magnetic devices. Some preliminary results have been reported already?*^ 

The paper is organized as follows. In Sec. II, we describe our system, the magnetic 
cantilever, and derive a set of equations describing magnetomechanical motion. In Sec. Ill, 
we solve these equations in the limit of small magnetization oscillations, finding a splitting 
of FMR spectra for the resonantly coupled mode. In Sec. IV, we propose a current-induced 
magnetic resonance technique in spin valves and demonstrate that magnetovibrational modes 
can be detected by this method. Finally, in Sec. V, we analyze large magnetization os- 
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cillations in the presence of magnetovibrational couplings and demonstrate magnetization 
reversal assisted by the mechanical degrees of freedom. 



II. SYSTEM 



We consider first a small dielectric cantilever with a single-domain ferromagnetic layer 
deposited on its far end (see Fig. [TJ). A constant external field H is applied along the z 
axis. In section 3, we also include an oscillating field H y along the y axis to pump and probe 
the system. The effective field H e fj felt by the magnetization consists of H as well as the 
crystal anisotropy and the demagnetizing fields. The strains are localized in the mechanical 
link of the cantilever between the ferromagnetic film at one end and the other end that is 
fixed. The lattice of the ferromagnet then oscillates without internal mechanical strains, 
but crystalline and form anisotropies couple the magnetic order parameter to the torsional 
mode of the cantilever. 

Our setup consists of two weakly interacting subsystems - the magnetic and the mechan- 
ical. It is then useful to first study the two subsystems separately. 



A. Magnetization motion in a single-domain ferromagnet 

The magnetization M of the ferromagnet precesses around an effective magnetic field 
H cff according to the Landau-Lifshitz-Gilbert equation^ 

TT a dM 

_ = - 7 MxH« ff + -Mx— , (1) 

where 7 denotes the gyromagnetic ratio. The phenomenological Gilbert constant is typically 
a < 0.01 for metallic and a < 0.00001 for insulating ferromagnets. The effective field is 
given by the functional derivative of the free energy of the system, that is of the form (in 
the lowest order in magnetizations in some specially chosen reference frame) 

E mg = X -D x Ml + \d v MI + h) z Ml - MH , (2) 

where D x , D y and D z describe the anisotropy of the magnetization along the Cartesian 
axes x, y and z, including demagnetizing effects and crystalline anisotropy. The associated 
effective field is given by H c g = H — D x M x y. — D y M y y — D z M z z. 
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B. Small magnetization oscillations and dependence of FMR broadening on shape 
and crystal anisotropics 



To first order, the deviations from the equilibrium magnetization in the z-direction lie in 
the x — y plane: M = m x x + m y z + M s z, where M s is the saturation magnetic moment. The 



magnetic susceptibilities Xy y (^) = ( m y/H y ) w and % 



(m y /H x ) describe the linear 



response of the magnetization m y to a (weak) rf magnetic field H y (H x ) at frequency u and 
can be found after linearizing the LLG equation in frequency space: 

7 2 M S (# + (D x - D Z )M S ) 



;i + a 2 )< + iau(2H + (D x + D y - 2D Z )M S 



(3) 



icujM s 



W > ~ ^ _ (! + a 2^2 m + muJ (2H + (D x + D y - 2D Z )M S ) ' 1 ] 

where u 2 m = (H + (D x - D Z )M S )(H + (D y - D Z )M S ) is the FMR resonant frequency. 
Note that the FMR damping has to be renormalized in the presence of anisotropies. The 
broadening of the FMR lineshape is 

(H /M s + (D x + D y )/2 - D z ) 
a ^(H /M s + D X - D z ) (H /M s + D y -D z ) 

Let us consider an ellipsoid with the semi-axes a and b = c. For small crystalline 
anisotropies the anisotropy factors are defined by shape anisotropies 



At 
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- \J m? — 1 ' 



4/T — 27rm, 



~ 7rm, 



(5) 



where m 
Fig.0 



a/c. The dependence of FMR broadening on the aspect ratio m is plotted in 



C. Nonlinear magnetization oscillations 

Without Gilbert damping the dynamics of the magnetization can be analyzed analyti- 
cally for different magnetic anisotropics^. We present here solutions for the "easy plane" 
anisotropy when D y = D z = and H = H z which is relevant for Section 4. Without the 
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Gilbert damping but including the anisotropy factor D x the LLG equations read 



^ = —jMyH 

< ^ = <yM x H + 1 D x M z M, 



X 



(6) 



^ = 1 D x M y M x 



V 



The x— component of the magnetization then obeys 



M x 



(H 2 - D x E m9 )j 2 M x - ^D\Ml 



(7) 



where the energy E mg = —H M Z + D x M 2 /2 is constant when there is no Gilbert damping. 

This equation describes a so-called "Duffing oscillator", one of the rare examples of a 
non-linear dynamic systems that can be solved analytically. The potential energy of the 
oscillator has two minima. As a result, one can observe effects like periodic motion centered 
at one of those minima that with increasing energy suddenly doubles its period. When a 
periodic external force is applied to such an oscillator, the system may carry out jumps 
between minima leading to stochastic motion in time. Duffing's equation 



where Z is the energy of the Duffing oscillator. By comparing Eqs. ()6I7|) to Eqs. (|8I9|) we 



obtain that Z = 7 2 (M 2 # 2 - E 2 ), k = ^ 2 (D x E mg - H 2 ) and A = 7 2 £> 2 /4. The energy 



minima of the Duffing oscillator therefore correspond to the energy maxima in our system. 
One can distinguish three different types of solutions. 

1. Z < 0. "Small amplitude" solutions when motion is at energies close to an energy 
minimum (in Fig. 01 it corresponds to two small circles and E mg > 1, the magnetization 
oscillates close to the perpendicular to the film direction). We express the solutions in terms 
of Jacobi elliptic functions dn: 



x - kx + 2\x 3 = 



(8) 



can be integrated at once to give the first integral: 



x 2 — kx 2 + Xx 4 = Z 



(9) 




(10) 



where v can be found from Z 



4A(1-d 2 /2) 2 ' 
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2. Z = corresponds to the separatrix (in Fig. 01 it is the longest possible trajectory in 
a shape of a bent "8"). Here the Jacobi functions degenerate to hyperbolic function ch: 

fk 1 

3. Z > (in Fig. EI it corresponds to bent elliptical trajectories and E mg < 1, the 
magnetization trajectories are squeezed in the perpendicular to the film direction). "Large 
amplitude" solutions via the Jacobi functions cn: 



x = ^w^) vcn [ V w^T) M 1 (12) 



k 2 v 2 (l-v 2 ) 



where Z - - A(2t)2 _ 1)2 

The analytical solution in our case can be written for the entire parameter space with 
the exception of the separatrix by a single formula as follows: 



M x 



^=2p~Zdn (t 7v /-p_/2, 1 - p + /p^ P ~~ P+ ^=2pZcos{t^/p~j2) (13) 



M — _ M x M _ D x Ml-2E mg (ia\ 



where p± = H$ - D x E mg ± H oy > ' H o ~ 2D x E mg + D 2 X M*. When D x > this leads to the 
trajectories and oscillation periods depicted in Fig. 03 The motion is periodic with time 
periods that can be expressed by elliptic integrals K as 



T x = 4V2K (p-/p+) I ( 7v ^T) P ~^ P+ 2n^2/p^; E mg < MH 



T 2 = 2V2K (1 - p + /p_) I ( 7 y=pT) lP ~ P ^ <<lP ~ l 7i^2/p~. E mg > MH 

The variation of the periodicity has important consequence for the coupling to the lattice in 
the regime of large magnetization oscillations as explained below. 

D. Cantilever oscillations and coupled magneto-mechanical equations 

Throughout this paper, we assume that the torsional oscillations of the cantilever are 
small. The torsional motion of the part of the cantilever that is not covered by the ferro- 
magnet can be found by applying the variational principle to the total elastic energy^ 

E el = \ ["cr'dy, (16) 



(15) 



o 
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where r = dip/dy and C is an elastic constant defined by the shape and material of the 
cantilever (C = ^pda 3 for a plate with thickness a much smaller than width d, a <C d, n 
is the Lame constant). T = Cr(y) is the torque flowing through the cantilever at point y. 
The integration is taken from the clamping point y = until the cantilever endpoint y = L. 
The equation of motion reads 

C%--^ + W% (17) 

where / = j(z 2 + x 2 )dzdx ~ ad 3 / 12 is the moment of inertia of the cross-section about its 
center of mass, p the mass density, and j3 is a phenomenological damping constant related 
to the quality factor Q at the resonance frequency oo e as Q = u e /(2/3) (at 1 GHz Q ~ 500). 4 
Note that u e can be also a higher harmonic resonance frequency in what follows. The 
oscillating solution has the form <p = sm(ky)(A 1 sm(ut) + A 2 cos(ut)), where k = (uj + i(3)/c 



is the wave number, c = c t 2h/d = C / (pi) and c t = \fp~Jp is the transverse velocity 
of sound. The free constants A\ and A 2 depend on the initial conditions. The boundary 
condition <p\ y =o = at the clamping point is already fulfilled, and the boundary condition 
at the end y = L is discussed in the following. 

Combining Eqs. ()2|16|) and taking into account the smallness of the magnet AL <Iwe 
can write the free energy of the cantilever coupled to the magnetization: 

F = y(-MH„ + -± [M x + MML)} 2 + ~y [M m - M X ^(L)] 2 + - JL-L) + - / r 2 dy. (18) 

J 

Eq. (|T8|) demonstrates that magneto-mechanical coupling is only possible when the factors 
D x or D z are non-zero (the anisotropy factor D y does not contribute to the coupling since 
our mechanical motion is rotationally invariant with respect to the axis y and one needs 
anisotropies D x and D z to break the invarance). For small y? = (p(L), where ip(y) is the 
torsion angle at position y of the cantilever, the effective field H cff = — is 

H efr = (D x M z (p - D x M x )x + D x M x (pz - (D z M xV + D z M z )z - D z M zV x - D y M y + H 

where D x , D y and D z describe the anisotropies (D x ~ 4tt for a thin film without crys- 
tal anisotropies). The equation of motion of the magnetization M in the presence of the 
mechanical degree of freedom then reads 

dM a f dM 



dt 1 M s \dt /aBt 



7 M x H cff + — M x — , (19) 



where the derivative (^r) cant = IF + T^ ( — M z x+M x z) is taken in the reference system of the 
cantilever since the magnetization damping is caused by interactions of the magnetization 
with the bulk of the cantilever. 

By applying the variational principle to Eq. (fT8f we obtain Eq. (fTTj) and the second 
boundary condition for its solutions. The magnetovibrational coupling can then be treated 
as a boundary condition to the mechanical problem, which is expressed as the torque Ct\ v= l 
exerted by the magnetization on the edge of the cantilever: 

CV| y=i = i P^+ T MxH J \ y , (20) 

This boundary condition is equivalent to the conservation law of the mechanical angular 
momentum written for the tip of the cantilever. Let us introduce an angular momentum 
V el (y) for a thin slice at point y G {0, L} (without magnetic overlayer), then the conservation 
law is <iV el (y) jdt = T (y) , where T(y) is the torque flowing into the slice. This equation is 
modified by the coupling to the magnet in a region y G {L, L + AL} (AL is the length of 
the cantilever covered by the magnetic layer) as: 

j f (v cl (L) + (-i)M (L) V^j = T (L) + T field . (21) 

where T field = V"M x H , V is the volume of the magnet and T(y)\ y = —Cr(y). When 
AL L, internal strains in the magnetic section may be disregarded and the magnetization 
torque can be treated boundary condition Eq. (l20j) . 
The coupling of Eqs. JTJ) and (jT7j) can be made explicit: 

H cff = (D x M z (p - D x M x )x + D x M x tpz - (D z M x <p + D z M z )z - D 2 M zV x - D y M y + H 

(22) 

E. Entirely ferromagnetic cantilever 

It is straightforward to generalize the free energy Eq. (JTSl) to the case when the whole 
cantilever is covered by a ferromagnetic film or when the whole cantilever is ferromagnetic. 
The free energy then has the form: 

F = V(-MH + ^ [M x + MM 2 + ^[M Z - MM 2 + _!L_L) + ± J r 2 dy (23) 
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where ip = jr J Q tp(y)dy is the average angle of torsion. In this case the system of the 
magnetomechanical equations becomes: 

f = - 7 MxH eff + |Mx(f) cant 
C 9 ^ = Plfg + f + 7M x H ) + 20pl% 
Cd(p/dx\ x=L = 

H cfr = (D x M z (p - D x M x )x + D x M x ipz - (D z M x ip + D z M z )z - D z M z ipx - D y M y + H 

(24) 

The coupling can be much more efficient in the latter as compared to the former system 
since the volume of the magnet is larger. The disadvantage is that the form anisotropy wants 
to align the magnetization along the cantilever making the subsystems uncoupled. Another 
disadvantage, the possibility of spin wave excitations, is discussed in the next section. 

III. SMALL MAGNETIZATION OSCILLATIONS 
A. Magnet at the tip of cantilever 

To first order, the deviations from the equilibrium magnetization in the ^-direction 
lie in the x — y plane: M = m x x + m y y + M s z, where M s is the saturation magnetic 
moment (see Fig. 1). For small ip, the effective field oscillates in the x — y plane: 
H cfr = (D x M B <p - D x M x )x + D x M x <pz - (D z M xV + D z M s )z - D z M s <px - D y M y + H , 
(D x ~ 47r for a thin film without crystal anisotropy). Note that the coupling does not rely 
on a strong magnetocrystalline field, since the surface forces of demagnetizing currents also 
provide a restoring torque. For a thin ferromagnetic layer the latter dominates and the 
magnetization does not precess, but oscillates like a pendulum in the y — z plane due to the 
oscillating field in the x-direction H cff = D x M s ipx. 

By using Eq. (f20|) with the substitution of ip in r, Cr\ y=L = Ckcos(kL)(Ai sm(ut) + 
A 2 cos(Ljt)) = Chip cot(fcL), we obtain in frequency space and to first order in the magneti- 
zation oscillations: 

V 

Chip cot(kL) = — {—ium y — 'yH m x ) (25) 
--Cip^(u 2 + 2iPu-ul) } (26) 
where in the second line the cot has been expanded close to the resonance frequency u e = 
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c7r(l/2 + s)/L ( s is integer, here we concentrate mainly on coupling to the first harmonic 
and s — 0). 

The magnetic susceptibilities Xyyi^) = ( m y/-^y) w an d Xyxl 1 ^) = (, m y/H-a) w describe the 
linear response of the magnetization m y to a (weak) rf magnetic field H y (H x ) at frequency 
uj. Generalization of Eqs (I3|4j) in the presence of magneto- vibrational coupling can be found 
after writing LLG equation in frequency space with use of Eq. (|25jl , and taking into account 
smallness of the Gilbert damping a: 

, s = 1 2 M 2 S {D X — D z + H Q (1 - gtan(fcL))/M„) 

X ™ [UJ) ij 2 -u; 2 m + 2za'u;u; m + [u 2 -H (H + (Dy-D z )M s )]gt a n(kL) [ ' 



, \ = iujM s 

X y A ) UJ 2_ u 2 n + 2ia'uu m + [uj 2 - H {H + (D y - D Z )M S )} ^tan(fcL) 1 1 

where g = M 2 V(D X - D z )c 2 / (CLu 2 e ) = (2/n) 2 (D x - D z )(L/a) 2 (V/V c )(M 2 / //) is the magne- 

tovibrational coupling constant (V and V c are the volumes of the magnet and the cantilever 

respectively) and the unperturbed resonance frequency uj 2 ^ = (H + (D x — D Z )M S )(H + 

(D y — D Z )M S ). Note that the coupling constant is defined by the material parameters in the 

term M 2 / fj,, by the geometry in the term (L/a) 2 (V/V c ) and by both the geometry and the 

material parameters in the term (D x — D z ). The imaginary part of Xyyd^) is proportional 

to the FMR absorption signal. 

When the cantilever with the magnet is subjected to the rf magnetic field, the mag- 

netomechanical torques induce a mechanical motion. In effect we have then constructed a 

nano-scale motor that transforms the magnetic field oscillations into mechanical motion with 

amplitudes that are given by the susceptibility x^ y (^) = (W-^y)o, • The latter follows from 

the LLG equation in frequency space using Eq. (|2*HI) . and taking into account the smallness 

of the Gilbert damping a: 

iujgtan(kL) 

Xipy ^ = uj 2 - u 2 m + 2iaW m + [uo 2 - H (H + (D y - D z )M a )} gtan(kL) ^ 
In the absence of the external field Ho the resonance frequencies in the vicinity of the 
resonance frequency u> e can be found after expanding cot(A;L) as in Eq. (j2*fij) : 



^1(2) 



1/21 1/2 



uj 2 e +oo 2 m + gu 2 e ± ((K + u m f + gul)((u B - u m ) 2 + guj 2 e )) . (30) 



When the external field H is smaller than the demagnetizing field, Eq. (f30"|) holds with 

uj^ = (H + (D x - D Z )M S )(H + (D y - D Z )M S ). This fact can be used to tune the FMR 
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frequency in order to match the elastic frequency. In this limit of strong shape anisotropy 
(for example a thin ferromagnetic film with D x ~ 47r) and when H G/(M s j3ul) <C 1, Eqs. 
(jZzEH) simplify to 

XyyM ~ u 2 _ ^ + 2iaW m + u 2 g tan(JfeL) ' 1 j 

. . zdjgtan(fcL) . , 

^yH « — 2 , 0I — 2 : TTTv 32 

^ or — ti^ + 2za'u;c<; m + org tan(fcL) 

Imaginary part of x yy (^) corresponding to the rf absorption is plotted in Fig. [HI We observe 

a typical anticrossing behavior between an optically active and non-active mode, with level 

repulsion and transfer of oscillator strength. The intrinsic damping of the mechanical system 

(for MEMS Q factors can reach 10 4 , which quickly deteriorate with decreasing size, however) 

imposes an extra damping on the magnetization dynamics, which close to the mode crossing 

may dominate the intrinsic damping due to a small Gilbert constant a. In case (3/u > a' 

the damping growth when we move from purely magnetic motion into purely mechanical 

motion along one of the lines in Fig. 0J In case (3/u < a' the damping diminishes along 

such line. 

The amplitude of mechanical oscillations along one of the spectrum lines in Fig. 0] is 
plotted in Fig. 6. The efficiency of the nanomotor is maximal at resonance. 



B. Ferromagnetic cantilever 

Eq. describes magnetomechanical dynamics for a ferromagnetic cantilever or a di- 
electric cantilever covered by a ferromagnet in the whole. The LLG equation is exactly the 
same as in case of previous subsections apart from the fact that the torsion angle should be 
averaged over the cantilever. Following Eq. (|24l) and for frequencies close to the mechanical 
resonance, the coupling can be written in Fourier space as follows 

-C(p—(u 2 + 2if3u - uj 2 e ) = ^-{-iujm y - ^H m x ) (33) 

This result is identical to the expansion in Eqs. I|25l2fij) after substituting V c /2 by V. In the 
expression for the coupling constant g in Eqs (I27I28B we have V/V c = 1/2. The coupling 
constant is thus increased in comparison to the case of a small magnet at the tip of the 
cantilever in which F/F c < 1. 
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The density of a metallic single crystal cantilever (Fe) is higher p ~ 8000 kg/ni and the 
Lame constant \i ~ lOOGPa^ is of the same order as for Si. Consequently, the metallic 
cantilever has to be smaller in order to have the same resonance frequency with the Si 
cantilever. Most importantly, we have to fulfill the condition H > D Z M S that forces the 
equilibrium magnetization direction along the z— axis (alignment with the y— axis would 
result in zero coupling). This condition can be easily fulfilled in our geometry by not too 
strong magnetic fields since the largest demagnetizing factor is D x . 

It follows from Eqs. (|27l28jl that a ferromagnetic cantilever is more suitable for obser- 
vating the magneto-mechanically coupled modes since the coupling strength is larger by the 
factor V C /2V compared to the system with a small magnet at the tip of the cantilever. Close 
to the mechanical resonance frequency both systems will behave identically. 

C. Observation 

The magnetovibrational coupling in our cantilever is observable by FMR and calorimetric 
techniques,— but the signal of nanoscale magnets is small. It might therefore be preferable 
to detect the resonance by the static deflection of the same cantilever due to an additional 
constant magnetic field along the described by Lohndorf et alM. In our approx- 

imation the field H T creates a torque 7M y if T x, whose modulation at the FMR conditions 
should be detectable. 

Since the vibrational frequencies of state-of-the-art artificial structures are relatively low, 
the use of soft ferromagnets (such as permalloy), is advantageous. The magnetic mode 
frequencies are then determined by shape anisotropies. The FMR frequency and the me- 
chanical resonance frequency should not differ by much more than Au ~ for a pro- 
nounced effect. A Si cantilever with axdxL=(lx5x 50) [im (C = 10~ 13 Nm 2 ) has 
a torsional resonance frequency of the order of u e = 10 MHz. Taking our ferromagnetic 
layer of dimensions al x d x AL = 50nm x 5/im x 5fim (thickness, width and length), then 
Wy^ ~ 100 KHz, meaning that we should tune the magnetic resonance to u e ± 100 KHz 
to observe the "polariton". The necessary rf field if x depends on the viscous dampings of 
mechanical and magnetization motion. At low frequencies additional sources of damping 
complicate measurements^ and the coherent motion of the magnetization can be hindered 
by domain formation. Coupling to higher resonance modes^ 3 - or structuring of the ferromag- 
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net may help to carry out measurements. 

An actual observation of the predicted splittings would give information about e.g. the 
magnetic moment of the film and the broadening would yield the quality factor of the 
elastic motion. From a technological point of view the tunable damping due to the magne- 
tovibrational coupling might be interesting for optimizing switching speeds. A ferromagnet 
effectively absorbs microwaves and turns them into a precessing magnetization, which via 
the magnetovibrational coupling can be transformed into a coherent mechanical motion. 
On the other hand, the ferromagnet may interact with the mechanical motion, to cause a 
magnetization precession, which in turn emits polarized microwaves. The emission in the 
coupled regime is more energy efficient in comparison with a fixed magnetic dipole emission 
in the case of small Gilbert constant but relatively low mechanical quality factor. This might 
be interesting for e.g. on-chip communication applications. 

In this section we have calculated the magnetic susceptibility of a system with magne- 
tovibrational coupling by magnetocrystalline fields or via surface forces of demagnetizing 
currents. A condition for effective energy transfer from an external rf magnetic field into 
mechanical motion and vice versa has been established. FMR spectra are predicted to split 
close to the resonance, and to strongly depend on the mechanical damping. The predicted 
effects should be observable with existing technology, but a further reduction of system size 
would strongly enhance them. 

IV. PROBING MAGNETIZATION DYNAMICS BY SPIN-TRANSFER 
TORQUES 

Microwave fields are a perfect tool for probing magnetization dynamics, e.g. by FMR. 
However, with shrinking size of the magnets such methods become increasingly difficult 
since the power absorbed by the magnet becomes very small. In this section, we propose 
a method to probe the magnetization dynamics in which the role of the microwave field is 
taken over by spin-transfer torques. By driving an AC current through F|N|F spin valve 
structures with fixed and free layer magnetizations that are canted with respect to each 
other, we can excite the free layer motion (Fig. E() (excitation of the magnetization motion 
by DC currents has already been realized^) . The magnetization dynamics should in turn 
influence the resistance of the device that can be measured. When the magnetic layer has 



13 



a mechanically active extension, as in Fig. we can also detect the magneto-vibrational 
mode. 



A. Electrical detection of FMR 

We consider here F1|N|F2 multilayer structures with perpendicular magnetizations mj 
and m 2 assuming m 2 fixed by anisotropics or exchange biasing. An AC current exerts an 
oscillating torque T\ on the magnetization 

ri = h l + R n /R 1 -(3P 1 /P 2 

Io 2e(l + R u /R 1 )(l + R n /R 2 )-{3 2 1 ) 

with p = cos6, AR m = 1/G 1(2)T + 1/G mi - 2R n , 4i? 1(2) _ = 1/G 1(2)T - 1/G mh 
Pi(2) = R\(2)- / R\(2) , and 2R^i = 1/G\+, + l/G^th where Gri(2)t an d ^1(2)1 are conduc- 
tances of the left (right) ferromagnet including the left (right) normal layer, G\+\ and CrLi 
are mixing conductances of the middle normal metal with adjacent ferromagnet interfaces. 
The resulting magnetization dynamics causes oscillations of the resistance of the multilayer 
structure according to*^ 

m = R + Rl + Bt - + f-> 3 ± (1 - ± *S-(* ± *u», (35) 

(R^l + Ri)(Rn + R2) — (3 R1R2 

In the following, we consider only small deviations of mi from the perpendicular to m 2 

direction. All the equations can then be rewritten keeping only the leading terms with 

respect to small (3. The torque T\ is proportional to the current Jo in this approximation: 

°2e(l + R n /R 1 )(l + R n /R 2 ) 1 *' { ' 

where we defined an effective rf field H x along the axis x that creates the same torque as 
the oscillating AC current. Thus the response function (rriy/Io)^ = Kx yx (u) with 

KM S V = A_ 1 + (37) 

2e(l + R n /R 1 )(l + R n /R 2 y { ' 

and the susceptibility can be found from Eq. (pjj. To the first order in the mechanical 

damping constant (3 we can write the dynamical impedance as: 

Z(«) = fl + fl, + fl 2 ______ . (38) 

|_~(„) I normalized by the resistance of the locked perpendicular magnetizations is plotted in 
Fig. El The parameters of the spin valve are the same with the symmetric setup of ref.— The 
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current amplitude Iq is chosen to correspond to magnetization oscillations of 15 degrees which 
appears to be readily experimentally accessible. In Fig. [HI the FMR resonance corresponds 
to the dip in the absolute value of the impedance that should be easily detectable. 

B. Electrical detection of magnetovibrational mode 

The method described in the previous chapter can also be used for detecting the mag- 
netovibrational modes. The ferromagnet Fl is extended forming a ferromagnetic cantilever 
(Fig. HJ). The magnetization m 2 can have two directions (bold and dashed lines in Fig. EJ) 
that are equally suitable for this purpose. In principle, the direction plotted by dashed line 
can cause a constant torsion of the cantilever when a constant current is sent through the 
spin valve since in this case the spin-transfer torque is transferred to the mechanical torque 
via the shape anisotropy. We concentrate here on the direction of m 2 indicated by the bold 
arrow which can be easier realized in experiment. When an AC current is sent through the 
system the oscillating torque causes oscillations of the magnetization mi as explained in the 
previous section. At resonance, we can observe the magnetovibrational mode as a splitting 
of the dip in the absolute value of the impedance plotted in Fig. E3 The width of the dips 
is a measure of the damping that is approximately half of the width in Fig. |H1 for chosen 
parameters, which is consistent with the results of Sec 3.1 . 

In a long cantilever with a finite exchange stiffness, the macrospin magnetization motion 
can be complicated by spin waves. The effective field for the LLG Eq. ^ in the presence 
of nonuniform magnetization reads 

2A 

H cfr = -D x M x x - DyMy-y - D z M z z + — y 2 M + H (39) 

where A is the exchange stiffness. The lowest-energy bulk spin wave mode is along the 
cantilever with the frequency that can be found from the LLG equation written in Fourier 
space: 

-ioJ l(j£k 2 + (D y - D z )M a + H ) 

-7(|^ 2 + (D x - D Z )M S + Ho) -too 

where M = (m x x + m y y)e*^' +kr ^ + M s z. The resonance frequency is u sw = 
lJ(jfk 2 + (D x - D Z )M S + H )(^-k 2 + (D y - D Z )M S + H ). We estimate the difference 



m x /M s 
m y /M s 



(40) 
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between frequencies of the macrospin mode with k = and the longest wavelength 



mode with k = ir/L. For a thin film 4nM s jj-(ir/L) 2 « 0.2 GHz, where we adopted 
A = 2 x 10~ n J/m, M s = 1.4 x 10 6 A/m and L = 1 Bin. Since we excite the mechanical 
motion monochromatically by the AC currents, it is sufficient then to have the mode split- 
ting larger than the broadening. This shows that in principle we can design the magnetic 
and mechanical subsystems to avoid bulk spin wave generation provided the mechanical 
damping as well as the Gilbert damping are relatively small. Note that by applying an 
antiferromagnetic layer on top of the cantilever we can strengthen exchange stiffness thus 
diminishing the possibility of spin waves even further. 

It is also possible to include coupling of the mechanical motion to spin waves but this 
regime will not be considered here. 

V. LARGE MAGNETIZATION CONES AND MAGNETIZATION REVERSAL IN 
THE PRESENCE OF COUPLING 

In this section, we extend the linearized magnetization motion to large angles relevant 
e.g. for the magnetization reversal. 

A. Resonant magnetization oscillations and reversal 

We consider resonant oscillations of the mechanical and magnetic degrees of freedom. 
We restrict ourself here to the case when only one anisotropy direction is present (e.g. 
D y = D z = 0) which is relevant for very thin ferromagnetic films with and easy plane 
anisotropy when D x ~ 4ti. 



We first address the maximal possible coupling strength of a system described by Eq. (JUJ). 
Consider the two subsystems oscillating at a common frequency uj. The total mechanical 
energy is then E me = plLuj 1 kp\ 1 where y? is the maximal angle of the torsional motion. By 
equipartition this energy should be of the order of the magnetic energy E mg = M S VH . 
The maximal angle would correspond to the mechanical motion induced by full transfer 




The coupling of Eqs. |H) and (jT7j) is 




(41) 



H eff = (D x M z ip - D x M x )x + D x M x tpz + IX 
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of the magnetic energy to the lattice. By equalizing those energies, we find an estimate 
for the maximal angle of torsion (p = ^M S VH / (pi Loo 2 ). The coupling between the 
subsystems can be measured by the distribution of an applied external torque (e.g. applied 
by a magnetic field) over the two subsystems. The total angular momentum flow into the 
magnetic subsystem by the effective magnetic field is (M S V '/ 'j)u , whereas that corresponding 
to the mechanical subsystem at the same frequency is (pILujip ) uj. Their ratio is ip co/ (jH ). 
The maximum angle (p derived above is therefore also a measure of the coupling between the 
magnetic and mechanical subsystems. This estimate is consistent with the coupling strength 
of polariton modes at resonance of g = iplu 2 /(■yHo) 2 rs (f^D x M 8 / Hq in Eqs. (|27|28j) (in this 
estimate we consider the case of not too strong external fields when u w 7 \ZHqL\Ms) . An 
estimate for a cantilever with p = 2330 kg/m 3 (Si) and d = 100 nm (uj ~1 GHz) leads to 
g = iplD x M s / H ~ (L/a) 2 (M 2 /p) ~ 10~ 3 . Increasing (L/a) and M s or decreasing Lame 
constant p is beneficial for the coupling. 

Magnetization reversal by a magnetic field in the coupling regime can be realized even 
without any damping by transferring magnetic energy into the mechanical system. Since 
we find that <p <C 1 for realistic parameters, the subsystems undergo many preces- 
sions/oscillations before the switching is completed. The switching is then associated with 
a slow time scale corresponding to the global motion governed by the coupling or a weak 
damping relative to a fast time scale characterized by the Larmor frequency. The equation 
of motion for the slow dynamics (the envelope functions) can be derived by averaging over 
the rapid oscillations. To this end, we substitute LLG Eq. (T[T3|) and first of Eqs. (JHj), 
linearized in the small parameters a, (3 and <p Q , into the equations for the mechanical and 
the magnetic energies: 

dt^ me = — 2/3 i^me + Ct\ x =l dt \x=L, 

f t E mg = -H M Z + D X M X M X . 
We focus in the following on the regime H <^ D X M S , which usually holds for thin films and 
not too strong fields, in which the magnetization motion is elliptical with long axis in the 
plane and small M x even for larger precession cones. Disregarding terms containing higher 
powers of M x and averaging over one period as indicated by (...) 

(^) + (2f3E me ) = -VD x (M z M xl p), 
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By adiabatic shaping of time-dependent magnetic fields H (t) we can keep the two sub- 
systems at resonance at all times (H (t) does not change much since we do not consider 
large angle cones that are very close to the antiparallel configuration). The slow dynamics 
(fi(L, t) ~ y4(t)e*^ +7r//2 )* and M x ~ W(t)e lu}t in time domain is then governed by the equation: 

A + (3A = -gu> m /M B (-l + , 
W + a'uW = uM s (-l + %g^)A } 

we introduced a frequency u m = 7 a/ D x M s H (t) that at £ = coincides with the frequency of 
fast oscillations u. Substitutions <fi(L,t) ~ A(t)e wt and ~ W / (t)e*^~ 7r//2 ^ corresponding 
to 7r/2-shifted harmonics are also solutions, and the initial conditions determine the linear 
combination of two envelope functions, i.e. the beating pattern of two hybridized polariton 
modes 9 . When initially all energy is stored in one degree of freedom, M x is n/2 shifted from 
<p(L,t) and A 2 (t) = 0. Eqs. describe a (damped) harmonic oscillator with frequency 
y/gu m uj <C uj when D x W/4 <C MH Q (since u m does not change much we replaced it by 
averaged Cu m ). Such oscillatory behavior persists for general angles (except for motion with 
very large angle cones close to the antiparallel configuration). This is illustrated by Fig. HH 
which shows a numerical simulation of Eq. (|2*2"1) for an undamped system excited at t = 
by a magnetic field H at an angle 27r/3 with the initial magnetization. The number of 
periods necessary to transfer all energy from one subsystem to the other is therefore given 
by ~ 1/ (A^fg). Eq. also shows that for damping constants a > (p /ir or (3/u > ^/g/n 
the beating is suppressed. 

Fig. HOI illustrates that the mechanical system absorbs energy from the magnetic sub- 
system and gives it back repeatedly in terms of violent oscillations that are modulated by 
an envelope function on the time scale derived above in Eq. as plotted by dotted line. 
When the envelope function vanishes the magnetization is reversed and the systems seems to 
be at rest. However, since the energy is not dissipated, the momentarily silence is deceptive, 
and the beating pattern repeats. An efficient coupling requires that the frequencies of the 
subsystems are close to each other at each configuration, which was achieved in the simu- 
lation by the adiabatic modulation of the magnetic field H according to Fig. 2. However, 
the reversal process is robust; an estimate from Eqs. (HHl) for the necessary proximity of the 
resonant frequencies of the mechanical and the magnetic subsystems is Au ~ v /^;. In that 
case, the above estimates still hold. 
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Let us now consider the magnetization reversal by an antiparallel magnetic field H (t) 
that undergo slow change of amplitude to keep the subsystems at resonance. This method 
of magnetization reversal by oscillating demagnetizing field due to mechanical oscillations is 
strongly analogous to the reversal by time-dependent magnetic fields.^ However, in contrast 
to the time-dependent field reversal, we vary the field H (t) but not the rf fields so avoiding 
complicated time dependence of rf fields «i We can calculate the dependence H (t) in advance 
for a specific sample, or we can create feedback circuit by connecting metallic contacts to 
the magnetic film, thus monitoring the dynamics. Since we assume zero temperature, the 
dynamics in Fig. [TT] is initiated by assuming a tiny angle of magnetization at t — 0. We 
wish to illustrate here that making use of the magnetoelastic coupling can accelerate the 
reversal importantly. We can suppress the backflow of mechanical energy for example by 
a sufficiently damped mechanical subsystem, as evident in Fig. El by comparing the two 
curves for (3 ~ and (3 ~ 0.02 with vanishing Gilbert damping, a = 0. Alternatively, we 
may detune the external magnetic field out of the resonance precisely after the first reversal, 
effectively rectifying the energy flow from the magnetic into the mechanical subsystem. We 
observe that even without any intrinsic damping (a = (3 = 0) the unwanted "ringing" can 
be strongly suppressed (dashed line in Fig. EJ. 

The experimental realization of such magnetization reversal will be a challenge since the 
cantilever has to preferably work at high frequencies ~ 1 GHz. In the resonant reversal, a 
significant coupling strength of g ~ 10~ 3 requires that one tenths of the cantilever volume 
is a ferromagnet. We investigated here the non-linear dynamics of coupled magnetic and 
mechanical fields for a cantilever with a ferromagnetic tip. Employing the new dissipation 
channels, we propose new strategies for fast magnetization reversal and suppressed "ringing". 
We can make use of the additional mechanical damping or shape the external magnetic field 
pulses, thus quickly channeling-off magnetic energy when damping is weak. 

B. Non-resonant magnetization oscillations and reversal 

We propose a non-resonant mechanical reversal scheme analogous to "precessional" 
switching^. The effective field H e ff (see Eq. (|22| ) has a component perpendicular to the 
plane of the film H x ~ Mtp. Under a sudden mechanical twist this component acts like 
a transverse magnetic pulse about which the precessing develops. Alternatively, we can 
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suddenly release a twist preinstalled on the cantilever. This could be achieved by an STM 
tip bonded to an edge of the cantilever at (L, d) and pulling it up slowly to the breaking 
point. The mechanical response should be fast, i.e. react on a time scale (7<£>z/M) _1 , but 
there are no resonance restrictions now. We integrate the equation of motion numerically 
for a strongly damped cantilever j3/u ~ 0.15 initially twisted by <p — 0.2 and suddenly 
released at t = 0. We reintroduce an easy axis anisotropy described by DM z z. Adopt 
D = 0.05, a = 0.01 and no external fields. Fig. [12] displays the desired reversal. The rather 
severe overshoot, as in the case of the precessional switching technique, can be minimized by 
carefully engineering the mechanical actuation to be closer to the optimum "ballistic" path 
between M z — ±1. 

Summarizing, we demonstrate here a precessional reversal scheme based on the mechan- 
ically generated out-of-plane demagnetizing field without applied magnetic fields. In this 
reversal scheme, the sharp control of the resonance condition is not required, however, the 
cantilever has to be fast determining the reversal time. 



VI. CONCLUSION 



This paper reports a detailed study of magnetomechanical effects in small magnetic can- 
tilevers, summarizing some old results*^ and making several generalizations. We proved 
that for small magnetization oscillations close to the mechanical resonance frequency, a fer- 
romagnetic cantilever behaves like the dielectric cantilever with a small magnet at the tip 
analyzed before. Such a ferromagnetic cantilever has an enhanced coupling of mechanical and 
magnetic degrees of freedom and thus is more suitable for observation of magnetomechanical 
effects. A metallic ferromagnetic cantilever can be integrated into magnetoelectronic circuits 
as shown in Sec. IV. Such an integrated device has the potential as a fast transducer of me- 
chanical motion as an alternative to previous designs that take advantage of magnetomotive 
forces. In our case, the magnetomotive forces are supplanted by spin-transfer torques. Our 
strategy; that is to say to avoid magnetic fields, is similar to what happens in the field of 
random access memories in which there is a considerable effort to replace magnetic fields 
by employing spin-trasfer torques. The technique presented in this paper strongly relies on 
resonant magnetovibrational coupling in which the magnetic torques can be effectively trans- 
formed into mechanical torques and vice versa. In order to see magnetomechanical torques 
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in experiement, one needs to be able to handle small structures on micro and nano scale, the 
magnets have to be small enough to form a single domain. We beleive that the integration 
of magnetoelectronics and magnetomechanics is possible and should lead to devices with 
new functionalities. 
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Figure 1: A magnetomechanical cantilever supporting magneto- vibrational modes. On a dielectric 
substrate (such as Si) a single-domain ferromagnetic film is deposited at the free end. 

Figure 2: Dependence of FMR broadening on aspect ratio m = a/c for ellipsoid with the semi-axes 
a and b = c (Hq = M s ). 



1 M. L. Roukes, Phys. World 14 (2), 25 (2001). 

2 K. C. Schwab and M. L. Roukes, Phys. Today, July 2005. 

3 J. A. Sidles, J. L. Garbini, K. J. Bruland, D. Rugar, O. Zuger, S. Hoen and C.S. Yannoni, Rev. 
Mod. Phys. 67, 249 (1995). 

4 X. M. H. Huang, C. A. Zorman, M. Mehregany, and M. Roukes, Nature (London) 421, 496 
(2003). 

5 C. Friedsam, A. K. Wehle, F. Kuhner and H. E. Gaub, J. Phys. Condens. Matter 15, S1709 
(2003). 

6 K. L. Ekinci, Y. T. Yang, M. L. Roukes, J. Appl. Phys. 95, 2682 (2004). 

7 D. Rugar, R. Budakian, H. J. Mamin and B. W. Chui, Nature 430, 329 (2004). 

8 I. Bargatin and M.L. Roukes, Phys. Rev. Lett. 91, 138302 (2003). 

9 A. A. Kovalev, G.E.W. Bauer and A. Brataas, Appl. Phys. Lett. 83, 1584 (2003). 

10 A. A. Kovalev, G.E.W. Bauer and A. Brataas, Phys. Rev. Lett. 94, 167201 (2005). 

11 B. Hillebrands and K. Ounadjela (Eds.), Spin Dynamics in Confined Magnetic Structures, 
(Springer, Berlin, 2003); S. Maekawa and T. Shinjo (Eds.), Applications of Magnetic Nanos- 
tructures, (Taylor and Francis, New York, 2002). 

12 M. Jamet, W. Wernsdorfer, C. Thirion, V. Dupuis, P. Melinon, A. Perez and D. Mailly, Phys. 
Rev. B 69 (2), 024401 (2004). 

Figure 3: The period T\ (E mg < 1) and T2 {E mg > 1) of the magnetization dynamics in Eqs. (fl"5|) 
in units 2ir/ (Ho + D x M)Hq \ as a function of energy in units HqM. The inset shows a plot 
of typical trajectories at different energies on the unit sphere (D X M = WHq 

)• 



22 



Figure 4: Dependence of the resonance frequency of the coupled motion on the FMR frequency w m 
of the uncoupled magnetization ( g ~ 0.001). 

Figure 5: The oscillator strength corresponding to each resonance in Fig. 4 in arbitrary units 
plotted by full lines and their widths plotted by dashed lines (g ~ 0.001, [3/uj > a'). 

13 Th. Gerrits, H.A.M. van den Berg, J. Hohlfeld, L. Bar, and Th. Rasing, Nature 418, 509 (2002). 

14 C. H. Back, D. Weller, J. Heidmann, D. Mauri, D. Guarisco, E. L. Garwin and H. C. Siegmann, 
Phys. Rev. Lett. 81, 3251 (1998). 

15 H. W. Schumacher, C. Chappert, P. Crozat, R. C. Sousa, P. P. Freitas, J. Miltat, J. Fassbender 
and B. Hillebrands, Phys. Rev. Lett. 90, 0117201 (2003). 

16 J.C. Slonczewski, J. Magn. Magn. Mater. 159, LI (1996); L. Berger, Phys. Rev. B 54, 9353 
(1996). 

17 M. Tsoi, A. G. M. Jansen, J. Bass, W.-C. Chiang, M. Seek, V. Tsoi and P. Wyder, Phys. Rev. 
Lett. 80, 4281 (1998). 

18 T. Valet, unpublished. 

19 A. Yamaguchi, T. Ono, S. Nasu, K. Miyake, K. Mibu and T. Shinjo, Phys. Rev. Lett. 92, 
077205 (2004); G. Tatara and H. Kohno Phys. Rev. Lett. 92, 086601 (2004); S.E. Barnes and 
S. Maekawa, |cond-mat/0311039 . 



20 J-E Wegrowe, D. Kelly, Y. Jaccard, Ph. Guittienne, J-Ph Ansermet, Europhys. Lett. 45, 626 
(1999). 

21 J.Z. Sun, J. Magn. Magn. Mater. 202, 157 (1999). 

22 E. B. Myers, D. C. Ralph, J. A. Katine, R. N. Louie, R. A. Buhrman, Science 285, 867(1999). 

23 S. I. Kiselev, J. C. Sankey, I. N. Krivorotov, N. C. Emley, R. J. Schoelkopf, R. A. Buhrman and 
D. C. Ralph, Nature 425 (6956), 380 (2003). 

24 A. V. Kimel, A. Kirilyuk, A. Tsvetkov, R. V. Pisarev and Th. Rasing, Nature 429, 850 (2004). 

25 Z. Z. Sun and X. R. Wang, preprint (2005). 

26 T.L. Gilbert, Phys. Rev. 100, 1243 (1955); IEEE Trans. Magn. 40, 3443 (2004). 

Figure 6: Amplitude of mechanical oscillations along one of the resonances in Fig. The maximal 
efficiency of the "nanomotor" occurs when both subsystems are at resonance. 



23 



Figure 7: Electro-magneto-mechanical device. 



Figure 8: Dependence of the normalized impedance on the AC current frequency {a! = 0.02). 

27 C. Serpico, I.D. Mayergoyz and G. Bertotti, J. Appl. Phys. 93, 6909 (2003). 

28 L.D. Landau and E.M. Lifshitz, Theory of Elasticity, 2nd ed. (Pergamon, Oxford, 1959). 

29 M. G. Simmons and H. Wang, Single crystal elastic constants and calculated aggregate prop- 
erties, MIT Press, 1975. J. D. Bass, Elasticity of Minerals, Glasses, and Melts, pp. 45-63, in 
Mineral Physics and crystallography: a handbook of physical constants, edited by T. J. Ahrens, 
AGU, Washington, DC., 1995. 

30 J. Moreland, M. Lohndorf, P. Kabos, and R. D. McMichael, Rev. Sci. Instrum. 71, 3099 (2000). 

31 M. Lohndorf, J. Moreland, and P. Kabos, Appl. Phys. Lett. 76, 1176 (2000). 

32 G. F. Cochran, R. W. Qiao, and B. Heinrich, Phys. Rev. B 39, 4399 (1989). 

33 U. Rabe, K. Janser, and W. Arnold, Rev. Sci. Instrum. 67, 3281 (1996). 

34 A. A. Kovalev, A. Brataas, and G. E. W. Bauer, Phys. Rev. B 66, 224424 (2002). 

35 T. Ono, H. Miyajima, K. Shigeto, K. Mibu, N. Hosoito, T. Shinjo, Science 284 (5413), 468 
(1999). 

36 M. Tsoi, R.E. Fontana, S.S.P. Parkin, Appl. Phys. Lett. 83, 2617 (2003). 

37 D. Atkinson, D. A. Allwood, G. Xiong, M. D. Cooke, C. C. Faulkner and R. P. Cowburn, Nature 
Materials 2, 85 (2003). 

38 I. N. Krivorotov, N. C. Emley, J. C. Sankey, S. I. Kiselev, D. C. Ralph, and R. A. Buhrman, 
Science 307, 228 (2005). 

39 A.N. Cleland, Foundations of Nanomechanics (Springer- Verlag, 2002). 

40 P.B. Visscher and X. Feng, Phys. Rev. B 65, 104412 (2002). 

41 S. Wolf, D. D. Awschalom, R. Buhrman, J. Daughton, S. von Molnar, M. Roukes, A. Chtchelka- 
nova, and D. Treger, Science 294, 1488 (2001). 

42 S. S. P. Parkin, Applications of Magnetic Nano structures (Taylor and Francis, New York, 2002), 
p. 237. 

Figure 9: Dependence of the normalized impedance on the AC current frequency (io m = w ei 
a' = 0.02, (3/uj = 0.002). 



24 



Figure 10: Time- dependent response of the magnetomechanical system to an external magnetic field 
switched on at t = at an angle 27r/3 with the initial magnetization, in the absence of dissipation. 
Plotted are the f(L) and z— components of the magnetization (D X M = IOHq). 

Figure 11: Time-dependent response of the magnetomechanical system to an external magnetic 
field switched on at t = (D X M = WHo, a = 0). The importance of mechanical damping can be 
seen by comparing the dotted line ((3 = 0) with the full line ((3 = 0.02). The dashed line illustrates 
the dynamics when the subsystems brought out of resonance by reducing the external magnetic 
field to half of its initial value at t = 6.5 ((3 = 0). 
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Figure 12: Magnetization switching without external field by an initially twisted cantilever that is 
released at t = 0. The inset shows the corresponding magnetization trajectory on the unit sphere. 
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